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Accelerator driven system (ADS) subcritical nuclear reactors are under development around the 
world. They will be intense sources of free, 30-55 MeV gA decay at rest These ADS reactor 
neutrinos can provide a robust test of the LSND anomaly and a precise measurement of the leptonic 
CP-violating phase 5, including sign{cos{S)). The first phase of many ADS programs includes 
the construction of a low energy, high intensity proton or deuteron accelerator, which can yield 
competitive bounds on sterile neutrinos. 


With their ability to enrich spent fuel and their 
inability to melt down, accelerator driven system 
(ADS) subcritical reactors are being designed and 
prototypes are being built around the world. ADS 
machines in general will be powered by proton 
beams, of energies of 1-2 GeV with currents of 10- 
20 mA, which strike high Z targets. Beam energies 
of 400 MeV or more lead to the production of pi- 
ons. The 7r“ are generally absorbed in the target, 
while 7r+ stop and decay at rest (DAR), producing 
which also stop and decay at rest, producing 
with a known energy spectrum peaked between 30 
and 55 MeV. Like the Vg from ordinary critical re¬ 
actors, these are free byproducts of the reaction, 
requiring no modihcation to the reactor design and 
so are at no cost to the neutrino physicist. Early 
phases of ADS use less energetic beams which can 
produce Vg via the decay at rest of spallation iso¬ 
topes (IsoDAR). In this letter we describe some of 
the neutrino physics which can be done with these 
ADS reactor neutrinos. 

For example, in January 2015 the Belgian govern¬ 
ment approved the MYRRHA installation at a price 
of 1.5 billion euros, whose construction will begin 
already in 2017 and finish in 2021. It is an ADS 
research reactor with a maximum proton energy of 
600 MeV and current of 4 mA, and so in particular 
it will already be a powerful source of ^DAR 

For concreteness we will consider two configura¬ 
tions which are currently planned by the Chinese 
ADS (C-ADS) collaboration. The first is a proto¬ 
type called CI-ADS which consists of a deuterium 
beam with 100 MeV/nucleon and a current of 15 
mA. The second, corresponding to C-ADS phase II, 
is an 800 MeV proton beam with a current of 10 
mA. As the beam energy in the first setup is below 
400 MeV, it cannot provide /iDAR and so will 


be used for a Vg program which will constrain the 
sterile neutrino mixing angle Ogg. The second will 
use to constrain the sterile neutrino mixing angle 
providing a robust test of the LSND anomaly, 
and to measure the leptonic CP-violating phase S. 

These experimental goals and their imple¬ 
mentations are inspired by the DAEJALUS 
project [I]. However there are several key differ¬ 
ences. DAEJALUS uses multiple accelerators with 
various baselines to the detector. To determine 
which event came from which source, it requires that 
only one source runs at a time. Therefore the max¬ 
imum instantaneous current at each DAEJALUS 
source needs to be quite high, in the case of the 
far source they envisage 30-50 mA. On the other 
hand, ADS projects consist of a single source and 
so the beam can be truly continuous and the same 
number of events/baseline can be achieved with a 
much lower maximum instantaneous current. Sec¬ 
ond, DAEJALUS proposes the construction of its 
own cyclotrons, in fact this is necessary since no 
country’s ADS program envisage the construction 
of multiple machines with the 10-20 km separation 
required by their proposal. Thus in a DAEJALUS 
program, neutrino physicists need to pay for the neu¬ 
trinos. Third, as a commercial reactor ADS has a 
highly redundant design. Even the C-ADS inject¬ 
ing accelerators are redundant. This redundancy 
will allow ADS to have very little downtime. Fi¬ 
nally, in the case of C-ADS, the location in southern 
China has the advantage that the horizontal com¬ 
ponent of the Earth’s magnetic field is 0.38 G, as 
compared with 0.17 G at the original DAEJALUS 
location in South Dakota and 0.13 G at the second 
DAEJALUS location in the Pyhasalmi mine. As a 
result the main backgrounds, which arise from low 
energy atmospheric neutrinos, are reduced by about 
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a factor of two in our setup. 

In the first experimental setup, which is similar 
to that proposed in Refs. m El, we use the 100 
MeV/nucleon deuterium beam to knock neutrons 
out of ®Be. These neutrons are moderated and mul¬ 
tiplied by heavy water moderator and then absorbed 
by a 99.99% isotopically pure ^Li, producing ®Li. 
The eventual /3 decay of the ®Li produces our Vf. 
with an energy spectrum peaked at 6 MeV but ex¬ 
tending up to 13 MeV. These Ve travel isotropically 
and are eventually detected via the inverse beta de¬ 
cay reaction 

Vg + p'^—>■ n + e'^ (1) 

on Hydrogen in liquid scintillator detectors placed at 
distances between 5 and 30 meters. Neutrino oscil¬ 
lations between the three known flavors are inappre¬ 
ciable at distances below 100 meters, thus a deficit of 
detected Vg suggests oscillations to sterile neutrinos. 

Backgrounds are notoriously problematic in such 
experiments, and so it is essential to place multi¬ 
ple detectors at multiple baselines. The baseline- 
dependence of the neutrino disappearance is well- 
known for a single flavor of sterile neutrinos and so 
can be used both to confirm that the anomaly is 
truly caused by sterile neutrinos and also to break 
the degeneracy between the sterile neutrino mix¬ 
ing angle and the sterile neutrino mass splitting. 
In the case of C-ADS, obtaining multiple detectors 
is quite simple. C-ADS will be located in China’s 
Guangdong province. In the same province, Daya 
Bay’s eight 20 ton liquid scintillator detectors will 
no longer be needed. We therefore suggest using as 
many of these eight Daya Bay detectors as possible, 
although in our simulations we have conservatively 
estimated that only two will be available. 

Not only would the presence of these neutrino de¬ 
tectors be useful for neutrino physics, but the ADS 
collaboration is also interested in nearby neutrino 
detectors as, with their reasonably high neutrino 
event rate, they provide real-time monitoring of the 
ADS target station and so yield additional informa¬ 
tion about the reactor as well as a safety measure 
against, for example, an overheating of the core. 

The beam-off backgrounds in this experiment are 
quite well understood, as they are identical to those 
in current reactor experiments. They are essentially 
all in some way created by cosmogenic muons and 


so can be considerably reduced with vetoes near 
these muon events. We have included the back¬ 
ground from decays of spallation ®Li, which is so 
long lived that it cannot be effectively removed by 
muon vetoes at this depth. We have assumed a 
depth of 100 meters, which has been considered 
for the target station. This gives us 3 background 
events/detector/day, consistent with that observed 
by the Daya Bay experiment with a similar overbur¬ 
den. Although we have run a series of GEANT 4 
simulations of the target, described in Ref. [1], we 
have not considered beam-on backgrounds in this 
analysis, a third Daya Bay detector would be use¬ 
ful to confirm that this approximation is reasonable. 
However, with sufficient shielding even the low base¬ 
line of 5 m is many times the neutron interaction 
length and so we expect few neutrons from the tar¬ 
get station to arrive at the detector, especially if the 
detector is placed far from the forward beam axis. 

We consider a run with 5 x 10^® ®Li decays, each 
yielding an Vg. In the absence of sterile neutrinos 
this yields 2 x lO'^ IBD events in a 20 ton target 
volume with 12% free protons at a distance of 30 
meters. We assume a 5% uncertainty on this event 
rate and a 1% error per 50 keV bin on the back¬ 
ground rate. As we have multiple detectors, the final 
precision depends only weakly upon this uncertainty 
unless AM^ < 0.3 eV^, which is well outside the re¬ 
gion preferred by global fits of sterile neutrino data. 

Note that, since the ®Li spectrum is independent 
of the beam and target station parameters, for any 
deuterium beam energy and current and any tar¬ 
get configuration, one may obtain 5 x 10^® such de¬ 
cays for an experiment of some duration. Therefore, 
given any beam and target parameters, one need 
only estimate this runtime. Neutron production 
rates from 62 MeV and 200 MeV deuteron beams 
on ®Be have been measured in Refs. EE] where 
they were found to yield roughly 0.12 and 0.6 neu- 
trons/deuteron respectively. As has been stressed in 
Ref. [3] , a deuterium beam is advantageous in that it 
directly brings neutrons to the target, which are lib¬ 
erated via proton stripping and deuterium breakup. 
These fast neutrons have energies peaked around 
half of the deuterium beam energy and so they must 
travel a long distance before slowing enough to have 
a high interaction cross section with ^Li. This places 
strong constraints on the target geometry. For this 
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FIG. 1: The region of 0ee — AM^ parameter space which 
may be excluded with 2a (black), Su (blue) and Scr (red) 
of confidence by an IsoDAR run with a near detector at 
5 m and a far detector at 15 m. 


purpose we have considered a target with multiple 
^Li sleeves surrounded by a graphite reflector. Our 
GEANT 4 simulations indicate that nonetheless only 
about 6% of the neutrons are captured by ^Li in our 
particular target geometry. This means that we re¬ 
quire 8 X 10^"* neutrons and so 1.4 x 10^® deuterons 
on target. With a 15 mA beam this corresponds to 
5 years of running. However an optimization in the 
target geometry as in Ref. [7], so as to increase the 
neutron capture rate, could easily reduce this run¬ 
time. 

We approximate the fractional energy resolution 
of the Daya Bay detectors, 5 years from now, to be 

We have parametrized the Ve disappearance using 
the parameters 0ee and AM^ defined by 

Pee = 1 — sin^(20ee)sin^ ( 1.27 


AM^L'^ 


E 


(3) 


where L is the baseline in meters and E is the neu¬ 
trino energy in MeV. We have determined the pa¬ 
rameter values which may be excluded at 2cr, 3cr 
and 5tT for all values of pairs of baselines between 
5 and 30 meters in a 5 year run, using the crude 
approximation that the detectors are pointlike. 

Fixing the near and far detector baselines to be 
5 m and 15 m, our exclusions are shown in Fig. 
Already in this 2-detector setup, in the 0.5 eV^-100 
eV^ range this exclusion regime is competitive with 
the most sensitive proposed sterile neutrino experi¬ 
ments. This proposal is complimentary to reactor- 
based experiments that will measure 6ee because the 
®Li decay spectrum extends about 50% higher than 
the relevant part of the reactor spectrum, and so one 


has access to the same values of LjE with values of 
L and E which are nonetheless distinct from those 
so far probed, providing an additional degree of ro¬ 
bustness. With three to eight of the Daya Bay de¬ 
tectors the reach and robustness would be improved 
considerably. 

In the rest of this letter we will consider C-ADS 
Phase B with a 10 mA beam of 800 MeV protons. 
These create neutrinos via 7r+ and DAR 

7r+/r+-I-e+-b -b (4) 

The spectrum of /iDAR neutrinos is well known. 
To determine the normalization of the neutrino sig¬ 
nal we again suggest using several near detectors, 
such as Daya Bay detectors, now at baselines rang¬ 
ing from 30 m to 100 m from the target. The results 
that we report below use two Daya Bay detectors, 
one at 50 m and one at 100 m. The near detectors 
can determine the normalization using elastic elec¬ 
tron neutrino scattering. 

In addition, the near detectors will again play the 
dual role of a sterile neutrino search and a real time 
monitor of the target. However, the sterile neutrino 
search will now be most sensitive in the —>■ Vf, 

appearance channel. Thus it is sensitive not to the 
sterile neutrino mixing angle 0ee but to defined 

by 

P^e = sin^(20^e)sin^ ^1.27 ^^ (5) 

and so it will test of the LSND anomaly [5]. Like 
LSND, as an appearance experiment its reach will 
extend to much lower values of the relevant mix¬ 
ing angle than a disappearance experiment. In fact, 
as the accelerator current is an order of magnitude 
greater than LSND and as there are multiple detec¬ 
tors, it will provide a much more robust test of this 
anomaly than LSND itself, ranging down to lower 
mass splittings and smaller mixing angles. 

The 2tT, 3cr and 5a bounds which may be achieved 
on and AM^ in a 5 year run are presented in 
Fig.i Note that the best fit LSND anomaly region 
is well covered, and so this will provide a robust test 
of the LSND anomaly. Elastic e — v scattering in 
the near detectors can be used to determine the 
flux with an uncertainty of about 5%, which is es¬ 
sentially irrelevant for such an appearance channel 
sterile neutrino search. The beam-off backgrounds 
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FIG. 2: The region of — AM^ parameter space which 
may be excluded with 2a (black), Su (blue) and Scr (red) 
of confidence by a 5 year /rDAR run with a near detector 
at 50 m and a far detector at 100 m. 


are dominated by atmospheric neutrinos which will 
be relevant at the far detectors but at the near de¬ 
tectors yield backgrounds which are again negligible. 
The beam-on backgrounds are more difficult to an¬ 
ticipate, but little is expected in this energy range. 
Following the estimates given by LSND [8] we fix the 
beam-on background by asserting that for each /i“'" 
DAR there are 5 x 10“^ DARs and that the nor¬ 
malization of this background is known with a preci¬ 
sion of 25%. We find that this beam-on background 
has little effect on our results. If there does happen 
to be an unexpected beam-on background, its con¬ 
tribution to the appearance rate at the near detec¬ 
tors would scale as 1/r^, in contrast with the sterile 
neutrino appearance signal, and so the two could 
in principle be disentangled. Of course in such a 
case it would be reasonable to bring more Daya Bay 
near detectors at different baselines to distinguish 
between the two signals as robustly as possible. 

So far the proposed physics program has been 
quite cheap, as the v have been supplied for free 
by ADS, albeit with a modified target station in 
the first setup, and the detectors have been recycled 
from Daya Bay. However perhaps the most exciting 
potential application of ADS reactor neutrinos is the 
measurement of the leptonic CP-violating phase 5. 
For this purpose, we will consider two 20 kton liq¬ 
uid scintillator detectors at baselines of 2.5-30 km to 
detect Ve created by oscillations. Combining this 
appearance result with —>■ Ve from NOi^A and 
T2K one obtains the leptonic CP-violating phase 5. 
Unlike the 20 kton detectors needed for mass hier¬ 
archy experiments, the measurement of 5 does not 
impose stringent restrictions upon the detector’s en¬ 
ergy resolution or the knowledge of its nonlinear en¬ 
ergy response. Thus, unlike JUNO [9], the detector 


can be cylindrical, the PMT coverage can be less 
than 50% and the scintillator purity does not need 
to exceed that in present day experiments. Also, as 
80% of the events are at energies above 30 MeV and 
so well above the spallation isotope backgrounds, the 
detectors do not need to be nearly as deep as JUNO. 

Such a setup has numerous advantages for mea¬ 
suring 5 nniiii]. Not only is the iZg spectrum well 
known, but since they are detected via IBD the cross 
section is also well known and the capture of the re¬ 
sulting neutron provides a double coincidence which 
can be used for background rejection. As a 17 oscil¬ 
lation experiment it maximizes synergy with beam 
experiments such as T2K and NO:/A, which are most 
efficient in neutrino mode. The V energies are higher 
than geoneutrino, reactor neutrino or even spalla¬ 
tion neutrino energies, and yet low enough so that 
atmospheric neutrino backgrounds are small. Un¬ 
like beam experiments the V energy is also below 
the muon mass, and so there is no background from 
charged current or interactions. The spec¬ 
trum is broad enough that such experiments can 
cleanly distinguish 5 from tt — 6. Finally the main 
background that appears in water Cherenkov ^DAR 
experiments mm, the invisible muons, are not 
present in a scintillator detector as even low energy 
muons, via ionization, emit scintillation light. 

We have considered a 12 year /rDAR run together 
with 12 years of NO^A appearance data, half in ly 
mode and half in V mode, defined as in Ref. m- 
The main background comes from atmospheric Ve 
and z7e interacting via IBD and also charged current 
quasielastic interactions. Roughly adapting the re¬ 
sults of the simulations in Ref. m, we determine 
the shape of this background and estimate the to¬ 
tal number of events to be about 50 in a 12 year 
run in southern China. In other locations the back¬ 
ground rate will be higher, as the strong horizontal 
component of the magnetic field in southern China 
leads to a smaller low energy neutrino flux. With 
two 20 kton detectors, one at 5 km and another at 
25 km, the precision with which <5 can be measured 
is given in Fig.[^ The precision, about 15°, is com¬ 
petitive with other future proposals such as DUNE 
and Hyper-Kamiokande but at a lower cost. Even 
with a single detector, as can be seen in Fig. a 
precision of about 25° can be expected. Needless 
to say, a single detector can be built first and the 
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FIG. 3: The la precision with which S can be measured, 
for each true value of S, by a 12 year /iDAR run with 
Daya Bay near detectors and 20 kton liquid scintillator 
detectors at 5 km and 25 km. 12 years of appearance 
data at NOi^A is also used. The dots are our numerical 
results, whereas the curve is an interpolation. 
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FIG. 4: The 5-averaged Icr precision with which S can 
be measured by ^DAR with a single 20 kton detector at 
the given baseline if the true value of 5 is 0° (black), 90° 
(red), 180° (blue) or 270° (green). Solid curves include 
NOr/A. 


second detector can be built later if and when the 
necessary funds present themselves. 
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